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We use entanglement to investigate tiie transition from vortex liquid phase to vortex lattice phase 
in weakly interacting rotating Bose-Einstein condensate (BEC). For the torus geometry, the ground 
state entanglement spectrum is analyzed to distinguish these two different phases. The low-lying 
part of ground state entanglement spectrum, as well as the behavior of its lowest level change clearly 
when the transition occurs. For the sphere geometry, the entanglement gap in the conformal limit 
(CL) is also studied. We also show that the decrease of entanglement between particles can be 
regarded as a signal of the transition. 
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I. INTRODUCTION 

Recently trapped rotating ultracold atoms have at- 
tracted considerable interest [i|, 0]. When the rotating 
frequency is moderate, the formation of vortex lattice 
is observed in experiments. However, with the increase 
of rotating frequency, it is predicted that vortex lat- 
tice will melt and the system will enter into a regime 
of strongly correlated many-body physics, namely vortex 
liquid phase [3, Q . The ground states in this liquid phase 
are the bosonic version of quantum Hall states, leading 
to the invalidity of mean field theory [5^9] . The parame- 
ter controlling this quantum phase transition from vortex 
lattice phase to vortex liquid phase is the average filling 
factor v = N/Nv, where N is the number of bosons and 
Nv is the number of vortices. Exact diagonalization on 
the torus geometry provides an evidence that this transi- 
tion occurs at the critical filling factor Vc ^ 6- For i^ > j^c, 
the ground state is a vortex lattice while for v < Vc, the 
ground state is a strongly correlated liquid [J]. 

The description of condensed matter phases using 
quantities such as entanglement entropy and fidelity 
borrowed from quantum information theory is another 
promising research field attracting great attention [lOJ. 
It has been found that when probing topologically or- 
dered states, topological entanglement entropy, which 
is the sub-leading term of the ground state entangle- 
ment entropy, can provide some information that cannot 
be obtained by conventional condensed matter methods 
pH [l2| . Quantum Hall states of two-dimensional elec- 
trons in magnetic field and two-dimensional bosons in 
rotating traps are very important topologically ordered 
states, whose topological entanglement entropies are cal- 
culated recently [13, JJ] . Comparing with generally used 
entanglement entropy which is only a single number, the 
entanglement spectrum (ES) contains more information 
and is found very recently to be useful in probing whether 
a quantum state has topological order [T^. In these two 
years, a series of papers focus on the ES of fermionic 



fractional quantum Hall states on the sphere geometry 
and torus geometry |15l - [l7| . the ES of topological insu- 
lators and superconductors [l^l, the ES of spi n model 
such as Heisenberg model and Kitaev model [19l42l| and 
the complete definition of entanglement gap in ES 1231 ■ 
However, the ES of rotating BEC still lacks enough study 
especially for large filling factors beyond quantum Hall 
regime and we will study it in this work. 

The ES is not only useful for systems of condensed 
matter physics. It is also important in quantum infor- 
mation theory. By majorization scheme of the ES, we 
can find deterministically whether a bipartite quantum 
state can be transformed to another by local quantum 
operations and classical communication while the entan- 
glement entropy is generally not sufficient [23| . 

In this work, we analyze the ground state ES of rotat- 
ing BEC to investigate the transition from vortex lattice 
phase to vortex liquid phase from an entanglement per- 
spective. Exact diagonalization method is used to ob- 
tain the ground state on the torus and sphere geometry. 
Then the ES is extracted from the ground state. Com- 
paring with Ref. J4|, where the excitation energy is used 
to identify the phase transition, we just utilize ground 
state properties, namely ES to indicate it. For the torus 
geometry, when bosons interact through contact inter- 
action, we find that with the increase of filling factor, 
the low-lying part of ground state ES as well as the be- 
havior of its lowest level undergo a qualitative change. 
The critical filling factor where this change happens is 
just Vc '^ (i- This phenomenon strongly indicates that 
a transition of ground state indeed occurs at Vc ~ 6. 
This is consistent with the conclusion of Ref. [J| . When 
we change the interaction to Coulomb potential, we find 
that the transition occurs at a little larger I'c- For the 
sphere geometry and contact interaction, the evolution of 
the entanglement gap in CL with filling factors is stud- 
ied. Finally, the entanglement between two bosons and 
other N — 2 bosons is also calculated. When ly > i/c, 
this entanglement decreases monotonically with i/, im- 



plying the transition from a strongly correlated regime 
to a mean field regime. 



II. MODEL AND METHOD 

Let us begin with a brief introduction of the definition 
of entanglement spectrum. Suppose we spatially divide 
a system into two parts A and B, then any pure state 
|*)yiB of the system can be expressed in Schmidt decom- 
position |^)ab = J2ie~^'^'^\4>i)A\ipi)B, where J2i<^~^' = 
1, £_i > 0, and {|(/'i)a} ({|<y5i)s}) is an orthonormal basis 
of subsystem A (B). The reduced density operator of A 

is pA = exp(-iJ) with H = J2i^^\4'i)A{4ii\- Considering 
in statistical physics a density operator can be related 
with Hamiltonian 7i through p = e~'^^/Tr(e~^^), We 
can regard H as a pseudo-Hamiltonian and define {^i} 
as entanglement spectrum. 

In this paper, we first study a rotating iV-boson sys- 
tem on the torus with periods a and b in the x and y 
directions. Each boson has a mass m. A consistent im- 
position of periodic boundary conditions requires that 
ab/£'^ = 2ttNv, where (. = \Jh/{2mu)) with a; the rotat- 
ing frequency and Ny is the number of vortex. When 
the interaction is weak enough, the lowest Landau level 
approximation is valid. The normalized single-particle 
lowest Landau level (LLL) wavefunction is 



( 1 V/2 ^ r./27r7 n6 

n^ — oo 

1 / . 2^J„2^2 



— (, + n5+— f 



where j = 0, 1, 2, ..., iVy — 1. Because V'j is centered along 
the line y = — ^^£^, the whole system can be divided into 
Nv localized orbits spatially. The interaction Hamilto- 
nian V{yi,T2) after standard second quantization pro- 
cedure becomes V = Efife~fc3fc4=o '^k^k2Hk^al^al^ak■,ak^, 
where a^ (aj,) annihilates (creates) a boson in the state 
i/'fe. On the torus, the matrix element Vtik-^Hki takes the 
form as 
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where 6' is the periodic Kronecker delta function with 
period Ny and V{qx,qy) is the Fourier transforma- 
tion of interaction V{r) [2J, ^^. For a Fock basis 

|no,ni, ...,njvv-i) satisfying I]jJo" % ^ ^' ^'^'^ *^*^^ 
find that there is a conserved quantity K — ^j^q' jrij 
(mod Ny), which can be called total momentum. The 
Lanczos algorithm is used to diagonalize the interaction 
matrix V to calculate the ground state energy and ground 



state wavefunction in each sector of K and then the low- 
est energy Eq and its eigenstate |$o) are selected. To 
extract the ES from \^q), we bipartition the system into 
blocks A and B which consist of Ia consecutive orbits 
and the remaining Ny — Ia orbits respectively. In this 
work, we focus on the case where Ia = Ny /2, namely A 
contains the orbits from to Ny/2 — 1 and B contains 
the orbits from Ny/2 to Ny — 1. ES is labeled by the 
particle number Na — X^^eA % ^^^ ^^^ total momentum 
Ka ~ J2j£AJ''^j (mod Ny) in block A. 

Similarly, we can consider a iV-bosons system on a 
sphere surface. Then under the LLL approximation, 
the contact interaction Hamiltonian can be written as 

V = Efcffc,fe3fe4=0^felfe2fe3fc44i42"fc3afc4 with 



Vk-i^k^k-iki — 0ki+k2Mi 



+k4 



1 (25+i)MntiC^f]^/' 



25 45 + 1 



^ki+k2 



where C", = („_,"^),„,, and Ny = 25. Uk («!) an- 
nihilates (creates) a boson with z-component angular 
momentum k. ES is labeled by the particle number 
Na = X^ieA "i ^^'^ ^^^ total z-component angular mo- 
mentum L^ = X^iga.^'^j ^^ block A. 



III. ENTANGLEMENT SPECTRUM 

Before we study the transition from vortex liquid to 
vortex lattice, we want to show some ES data of typi- 
cal bosonic quantum Hall states. We use two-body con- 
tact interaction J2i<j H^i ~ ^j) ^^'^ three-body contact 
interaction J2i<j<k^(''^i ~ i"j)^(''i ~ ^k) to obtain the 
exact Laughlin state at i/ = 1/2 and Pfaffian state at 
ly — 1 respectively. Throughout our calculation we take 
h/{muj) = 1, leading to ^^ — 1/2. When plotting the 
ES data, we choose the sector of Na where the lowest 
ES level is found. Then in this Na sector, the Ka sec- 
tor where the lowest ES level locates is put in the center 
of the spectrum. AKa is defined as the shift of mo- 
mentum from the center of the spectrum. For example, 
for Ny = 6, the possible values of Ka are 0,1,2,3,4 and 
5. If we find the lowest ES level locates in Ka = 2 
sector, then we plot ES data according to the order 
Ka = 5,0,1,2,3,4,5 from left to right to put Ka = 2 
sector in the center of the spectrum, corresponding to 
Ai^A = -3,-2,-l,0,l,2,3. 

The shapes of the ES of both Laughlin and Pfaffian 
states can be explained by considering their thin torus 
limit (a — >■ 0). We give some examples of such an anal- 
ysis below. For the Laughlin state with N — 2, Ny — 4 
in K — sector, the thin torus state is |0101), so the re- 
duced density operator of part A is pA — |01) (01 1 leading 
to only one level in the bottom of ES in Na — 1 sector. 
When we enlarge a, more ES levels come down and a V- 
shape appears [FiglIJa)]. For fermion Laughlin state, this 
V-shape can be explained by two chiral edge theory for a 
properly small a depending on the system size !17'|. For 



the PfafBan state with N = Ny = 4 in X = sector, the 
thin torus states are (|0202) ± |2020))/ V2, so the reduced 
density operator of part A is pA^ (|02> (02| + 120) (20|)/2, 
leading to two equal ES levels in Na = 2 sector with 
Ka = and 2 respectively. So the ES has a flat 
bottom with three degenerate levels in AKa = 0,±2. 
When we enlarge a, more levels come down but the 
flat bottom with three degenerate levels is still there 
[FiglUb)]. However, this flat bottom can disappear if we 
choose another system size. For the Pfaffian state with 
N = Ny = 6 in K — sector, the thin torus states are 
(1020202) ± 1 202020) )/V2, so the reduced density opera- 
tor of part A is PA = (|020)(020| + |202)(202|)/2, leading 
to two equal ES levels again. But they are in different 
Na sectors {Na = 2 and Na = 4 respectively). So, when 
plotting the ES in one Na sector, the flat bottom dis- 
appears and V-shape revives [Fig[TJd)]. At last, for the 
Pfaffian state with N = Ny = A in K = 2 sector, the 
thin torus state is |lf f f). So it is easy to see there is only 
one ES level in the bottom in Na — 2 sector, meaning a 
V-shape will appear when a increases [Fig[Hc,e)]. 

In Ref . '4s| , a trial wavefunction for vortex liquid phase 
at i/ = k/2 with k an integer is proposed: 



N/k 



^fe = s 



[ n (^^ 

i<jeAi 



N/k 

. n (^' 

l<rn^Ak 



where we express it in disk geometry and the Gaussian 
exponential factor is omitted. In every set Ai, there are 
N/k bosons and S means the symmetrization over all 
possible partition of the N bosons into k sets. When 
k = 1(2), this state is Laughlin (Pfaffian) state. We 
expect that the ES of incompressible quantum Hall states 
at other filling factors on the torus have similar shapes 
with Laughlin and Pfaffian states. Therefore if we see a 
V-shape or a flat bottom with three degenerate levels in 
the ES of one state, we can conjecture that this state has 
a close relationship with quantum Hall state. 

Now we consider a system on the torus with the aspect 
ratio a/b commensurate with a triangular lattice, which 
can help us to stabilize the vortex lattice. We suppose 
the interactions are contact interaction V{r) — 6{r) and 
Coulomb interaction V{r) = l/|r|, for which V{qx,qy) oc 
1 and V{qx, Qy) oc {q^ + q^^^l"^ respectively. Because we 
will consider the ground state at large filling factors, the 
number of vortex cannot be too large. First we study the 
cases at integer filling factors. When z/ is a small integer 
(e.g., V <\^ for Ny — 4), the ground state is unique and 
for our choice of aspect ratio it is in iiT = sector. When 
!/ is a large integer (deeply in vortex lattice phase, e.g., 
i^ > 10 for Ny — 4), the ground state energies in differ- 
ent K sectors become nearly the same. After checking we 
find all the ES of these nearly-degenerate ground states 
have qualitatively the same structure, so we still choose 
the ground state in K = ^ sector to extract the ES. In 
FigsH] and El we plot the ground state ES for Ny = 4 
and Ny = 6 at integer filling factors v — 1, 2, ..., 11, 12 
respectively. For Ny = 6, let us first consider contact 
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FIG. 1: The ES of typical quantum Hall states (the aspect 
ratio is 1). The Laughlin state is two-fold degenerate on the 
torus: one is in sector K = and the other is in sector 
K — Nv/2. The Pfaffian state is three-fold degenerate on the 
torus: two are in sector K — and one is in sector K = Nv/2. 
(a) The ES of Laughlin state at «/ = 1/2 for N = 8, Ny = 16 
in sector K — 0. (b) The ES of Pfaffian state aX u = 1 for 
Af = 8, iVv = 8 in sector K = Q. (c) The ES of Pfaffian state 
at !/ = 1 for Af = 8, A^v = 8 in sector K = A. (d) The ES of 
Pfaflian state at i/ = 1 for iV = 10, Ny = 10 in sector K = Q. 
(e) The ES of Pfaflian state &% v ^ I ior N = 10, Ny = 10 in 
sector K — 5. 



interaction. One can see that the qualitative behavior 
of the low-lying part of ES clearly changes after i/ = 6, 
distinguishing the vortex lattice phase and vortex liquid 
phase. For j^ < 6, the ground state ES has a structure of 
V-shape, while for v >7 , this V-shape structure vanishes 
and the ground levels of ES in every IS.Ka sector have 
a tendency to become degenerate, implying a transition 
from quantum Hall region. If we change the interaction 
to Coulomb interaction, the structure of the ground state 
ES changes after a little larger filling factor v — 'd. Simi- 
lar phenomena occur for Ny = 4. The qualitative behav- 
ior of the low-lying part of the ES clearly changes after 
V — 1 and V — 'd for contact interaction and Coulomb 
interaction respectively. Before the transition, the shape 
of the ES is either a V-shape (at even filling factors) or 
a fiat bottom with three degenerate levels (at odd filling 
factors). After the transition, the ground levels of ES in 
every AKa sector have a tendency to become degenerate. 
When the filling factor is a fraction -p/q where g ^ 1 
and p and q are coprime, the ground state is at least q- 
fold degenerate. For our choice of aspect ratio, there is 
always one of these degenerate ground states in K = Q 
sector. After calculating the ES of all degenerate ground 
states (e.g., for g = 2), we find that just like in the case of 
integer filling factors, with the increase of v^ the ground 
levels of the ES in every AKa sector have a tendency to 
become degenerate. 

In Ref. |J| , a ground energy gap of the Hamiltonian is 
defined which can indicate the transition from vortex liq- 
uid to vortex lattice (FigH]) [FigUJ^a) can also be found in 
Ref. [J]]. One can see for Ny ~ 6 the ground energy gap 
vanishes after i/ = 6 and v — % for contact and Coulomb 
interaction respectively, which is consistent with the con- 
clusion obtained by ground state ES. 
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FIG. 2: The ground state ES for Nv — 4 at integer filling factor i/ — 1,2, ..., 11, 12 from left to right. The subsystem A consists 
of orbits and 1 and the subsystem B consists of orbits 2 and 3. The Na sector is chosen as the one where the lowest ES 
level (corresponding to the largest eigenvalue of pa) locates. The aspect ratio is chosen as V3/2. The interaction can be either 
contact interaction (the first row) or Coulomb interaction (the second row). 
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FIG. 3: The ground state ES for Nv = 6 at integer filling factor i^ = 1, 2, ..., 11, 12 from left to right. The subsystem A consists 
of orbits 0, 1 and 2 and the subsystem B consists of orbits 3, 4 and 5. The Na sector is chosen as the one where the lowest ES 
level (corresponding to the largest eigenvalue of pa) locates. The aspect ratio is chosen as l/\/3. The interaction can be either 
contact interaction (the first row) or Coulomb interaction (the second row). 



Considering the critical point of a quantum phase tran- 
sition can be indicated by the singular behavior of the 
ground energy (or its derivative) of the physical Hamilto- 
nian and the ES corresponds to the spectrum of a pseudo- 
Hamiltonian H, we hope to see the ground energy of 
ES can indicate the transition from vortex liquid to vor- 



tex lattice. We investigate the evolution of the lowest 
level -Ecn of ground state ES with filling factor ly. Here 
we take a system containing six vortices as an example. 
For contact interaction, from FiglSJa) we can see when 
J/ < 5, Een increases slowly and oscillates obviously with 
ly. When j^ > 7, E'en also increases very slowly with i/ 
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unnormalization, the ground state is expressed as l"^) ex 
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FIG. 4: The ground energy gap A(A'^) of the Hamiltonian as 
a function of N for (a) contact interaction and (b) Coulomb 
interaction. A(A'') reflects the discontinuity of chemical po- 
tential and is defined as A(iV) = N\E[N + 1)/(A'' + 1) + 
E{N - \)I(N - 1) - 2E{N)IN\ where £;(A^) is the ground 
energy of the system consisting of TV particles. The vanish 
of A(A'^) indicates the transition from vortex liquid to com- 
pressible vortex lattice. Here Ny — 6 and the aspect ratio is 
l/vo. The unit of the energy is huj. 




FIG. 5: (color online) The lowest level Ecn of the ground 
state ES as a function of v for (a) contact interaction and (b) 
Coulomb interaction. When the ground state is degenerate 
(or nearly-degenerate), we choose the ground state in i^ = 
sector to extract the ES data. After checking, we find this 
choice will not influence our results qualitatively. The integer 
fllling factors are emphasized by red dots. Here Nv = 6 and 
the aspect ratio is l/y/S. 



E„o,"i,...,n2s('^"o,ni,...,"2slli=0-''^/il'^0,"-l,---,"-2S; 

where \n(),ni, ...,n2s) is the new basis state that does 
not contain A/j. The ES in CL is extracted from |^). For 
an incompressible quantum Hall state, the ES in CL can 
be divided into two parts, one of which is the low-lying 
conformal field part and the other is the generic part 
caused by realistic interaction. There is an entanglement 
gap between two parts [H, [22] • If a state undergoes a 
transition from a quantum Hall state to a compressible 
state, this entanglement gap closes [16l.[23|. 

Here we study the evolution of ground state ES in CL 
with the increase of filling factor on the sphere geometry 
(Figini). When the filling factor is small {v — 1), there 
is a clear entanglement gap separating the low-lying con- 
formal field theory part from the generic part. However, 
when the filling factor becomes large (i/ = 3/2, 2, 5/2, 3), 
the gap is unclear and hard to identify. Therefore there 
are not many hints for us to say wether the ground state is 
an incompressible quantum Hall state or not for v > 3/2 
in thermodynamic limit. Our conclusion obtained from 
the entanglement gap is consistent with the^one obtained 
from the energy gap of the Hamiltonian [26|] , where the 
authors found for v > 3/2 the energy gap is very small 
showing no clear signs of convergence towards thermo- 
dynamic limit and the two-particle correlation function 
shows very strong oscillations and no hint of incompress- 
ibility. The possible reason is that the quantum Hall 
states at larger filling factors have very large correla- 
tion lengths which exceed our largest finite-size system. 
Therefore, the ES, like the spectrum of the Hamiltonian, 
can detect the quantum Hall states at small filling fac- 
tors on the sphere. However, probably it is helpful to ob- 
serve a possible transition from an incompressible quan- 
tum Hall state to a compressible state only in much larger 
systems. 



but stops clear oscillation. A transition regime can be 
found at 5 < i^ < 7, in which E^n goes up steeply with 
the increase of v, showing a singular behavior like that in 
first-order quantum phase transition. Similar phenomena 
happen for Coulomb interaction (Figl5{b)). 

We wonder wether the transition from a bosonic 
quantum Hall state to a compressible state can also be 
refiected in other geometry. In Ref. [2^, the ES in the 
conformal limit (CL) on the sphere geometry is defined. 
The spirit of CL is to unnormalize the ground state 
obtained from the diagonalization of the Hamiltonian by 
expressing the ground state in a new set of basis states 
which do not contain the normalization factor of single- 
particle LLL wave function, because this factor contains 
the information about the extent of the Landau orbits 
in space and depends on the magnetic length of the 
problem. On the sphere, the normalization factor of j-th 
single-particle LLL is J\fj oc ^y{2S)\/[jl{2S - j)\],j = 
0,1,..., 25. If the ground state is |*) = 



Er. 



.'4>-n 



;|no,ni,...,n2s), 



then after 



IV. PARTICLE ENTANGLEMENT 

We know that in vortex lattice phase, the Gross- 
Pitaevskii mean field theory is essentially correct. The 
many-body ground wavefunction in mean field theory is 
just a product of wavefunctions of every particle, mean- 
ing that in mean field theory there is no ground state 
entanglement between particles. But in vortex liquid 
phase, the ground state is strongly correlated. So it's 
natural to ask whether particle entanglement can be an 
indicator of the transition from vortex liquid to vortex 
lattice. Here we consider the ground state particle en- 
tanglement between two bosons and other N — 2 bosons, 
defined as 52 = — Tr(p2lnp2)- P2 is the reduced den- 
sity operator of two bosons, whose matrix element is 
{P2)tj,ki = jv(j^_i) (aI"?"j°»)GS, where (•)gs is ground 
state average. Considering the case Ny = 6 on the torus, 
one can see that 52 decreases obviously after v — 6 for 
contact interaction and after v = 8 for Coulomb interac- 
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FIG. 6: (color online) The ground state ES in CL on the sphere geometry for contact interaction with filling factors v = 
1, 3/2, 2, 5/2 and 3 from left to right. Here A'^ = (Ny + l^v. Subsystem A consists of orbits from to Ny (2 — 1 and subsystem 
B consists of orbits from Ny /2 to A'^v'. Unlike what we do on the torus geometry, we choose the sector Na ~ N/2 for even A'^ 
[Na = {N - l)/2 for odd iV] to plot ES when considering sphere geometry, (a) z/ = 1, iV = 14, Ny = 12. (b) u = 3/2, iV = 18, 
Ny = 10. (c) u = 2, N = 20, Ny = 8. (d) u = 5/2, iV = 25, Ny =8. (e) z^ = 3, iV = 30, Ny = 8. The entanglement gap is 
clear at ;/ = 1 (indicated by the red lines), but hard to identify for larger v. 
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FIG. 7: (color online) The ground state particle entanglement 
52 on the torus for both contact interaction (black cubic) and 
Coulomb interaction (red circle) . Here Ny = 6 and the aspect 
ratio is l/-v/3. 



tion respectively (Fig|7]), giving a signal of the transition 
from vortex liquid to vortex lattice, though not very ex- 
actly. 



V. SUMMARY 

In summary, we investigate the quantum phase tran- 
sition from vortex liquid phase to vortex lattice phase 
in rotating BEC from an entanglement perspective. The 
ground state entanglement spectrum is a useful tool in 
detecting this transition. On the torus geometry its low- 
lying part has different structures in vortex liquid phase 
and vortex lattice phase. Its lowest level also behaves 
differently in these two phases. Therefore ground state 
ES distinguishes these two phases very well. Our results 
are consistent with the conclusion that for bosons inter- 
acting through contact potential, a quantum phase tran- 



sition happens at i/c ^ 6. When considering Coulomb 
interaction, i/c will become a little larger. On the sphere 
geometry, at i/ > 3/2, the entanglement gap of ES in CL 
is unclear, showing a need to investigate much larger sys- 
tems if one wants to observe the melt of vortex lattice. 
Finally, we show that particle entanglement can also be 
a signal to indicate the transition. 

In this work, we choose the aspect ratio as the one 
commensurate with triangular lattice. If we change the 
aspect ratio and the interaction, we can study the ES of 
other possible ground states, such as smectic state and 
vortex lattice state with other symmetry. Moreover, it's 
also interesting to investigate the ES of the simplest non- 
Abelian bosonic Moore- Read state at z^ = 1. In a very 
recent paper [23| , a new kind of ES is defined where spa- 
tial bipartition of the system is replaced by particle bi- 
partition. We hope this new ES can also indicate the 
transition from vortex liquid to vortex lattice. 

Current experiments are deep in the regime of vortex 
lattice, so it is still a big challenge to create strongly- 
correlated states in vortex liquid phase. Formation of 
these new states requires faster rotation and smaller par- 
ticle number than most current experiments. However, 
some method to solve this problem has been proposed 
[28J . We hope the transition from vortex lattice to vor- 
tex liquid may be detected experimentally in the future. 
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